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Abstract 

We propose a simple model in which the cosmological dark matter consists of par- 
ticles whose mass increases with the scale factor of the universe. The particle mass is 
generated by the expectation value of a scalar field which does not have a stable vacuum 
state, but which is effectively stabilized by the rest energy of the ambient particles. 
As the universe expands, the density of particles decreases, leading to an increase in 
the vacuum expectation value of the scalar (and hence the mass of the particle). The 
energy density of the coupled system of variable-mass particles ("vamps") redshifts 
more slowly than that of ordinary matter. Consequently, the age of the universe is 
larger than in conventional scenarios. 



1 Based on a talk by SMC at COSMO-97, International Workshop on Particle Physics and the Early 
Universe, 15-19 September 1997, Ambleside, Lake District, England. 
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1 Introduction 



The Big Bang model has proven extraordinarily successful as a framework for interpreting 
the structure and evolution of the universe on large scales. Within that framework, the cold 
dark matter scenario (featuring massive particles which bring the density of the universe 
to its critical value, and a scale-free spectrum of Gaussian density perturbations) has pro- 
vided an elegant theory of structure formation, which unfortunately seems to fall short of 
perfect agreement with observation. Although the precise extent to which CDM disagrees 
with observation is arguable, there are two important areas in which the discrepancies are 
particularly troubling: predicting an age for the universe which is larger than the ages of 
the oldest globular clusters, and matching the COBE-normalized power spectrum of den- 
sity fluctuations as measured by microwave background anisotropy experiments and direct 
studies of large-scale structure. 

One way in which the simple CDM scenario may be modified, affecting the age of the 
universe as well as the evolution of density fluctuations, is to imagine that the closure density 
is provided by something different than (or in addition to) nonrelativistic particles. In a flat 
Robertson- Walker universe with metric 

d s 2 = -dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ) (1) 

and energy-momentum tensor 

T^v = diag(-p,p,p,p) , (2) 
the Friedmann equations imply that the time derivative of the scale factor a satisfies 

•2 8nG 2 / Q \ 

a = —^-a p . (3) 

The evolution of a is therefore dependent on how the energy density p scales with a; if 
p oc a~ n , we have 

- = -—(2-n)p. (4) 
a 6 

Hence, the more slowly the energy density decreases as the universe expands, the more slowly 
the expansion will decelerate, implying a correspondingly older universe for any given value 
of the expansion rate today — for a flat universe dominated by such an energy density, the 
age is t = -Hq , where H = a/ a is the Hubble parameter and the subscript refers to 
the present time. (Eq. ([|) can also be derived by positing an equation of state p = wp and 
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using energy conservation; the two parameterizations are related by n = 3(1 + w).) The 
energy density in a species of ordinary "matter" (a nonrelativistic particle species X) can 
be written px = mxnx, where mx is the mass of the particle and nx is its number density. 
The energy density of a matter-dominated universe is therefore proportional to a -3 , as the 
mass stays constant while the number density is proportional to the volume; the age of such 
a universe is to = ^H^ 1 . 

Although there is some controversy over the value of the Hubble constant, most recent 
determinations are consistent with a value H = 70 ± 10 km/sec/Mpc, or Hq 1 = (14 ± 
2) x 10 9 yr [I]]. The upper limit on the age of a matter-dominated flat universe is therefore 
to(MD) < 11 x 10 9 yr. Meanwhile, calculations of the ages of globular clusters imply an age 
tec ~ 15 x 10 9 yr, with a lower bound of tec > 12 x 10 9 yr |2|. The apparent discrepancy 
between these values may be resolved by a revision in distance determinations to globular 
clusters, as suggested by recent measurements by the Hipparcos satellite ||; while this would 
be the simplest solution, further work is necessary to accept it with confidence. 

Alternative resolutions are provided by models in which the density parameter Q — 1, 
and some or all of the unseen energy density resides in a component which redshifts more 
slowly than nonrelativistic matter. The most popular such alternative is the introduction 
of a cosmological constant A, for which p\ = constant. Such models have some attractive 
features, but are also plagued with both theoretical and observational disadvantages (|, |J. 
A popular variation on this theme is to invoke a slowly-rolling scalar field, or equivalently 
a cosmological constant whose value varies with time, or simply an unspecified smooth 
component ||. More speculative possibilities include a network of cosmic strings or 
stable textures |§. We will not enumerate the good and bad qualities of each of these 
scenarios, noting only that none are sufficiently compelling to discourage the exploration of 
still further models. 

In this paper we propose a simple model in which the dark matter consists of particles 
ip whose rest mass increases with time. This is achieved by having the rest mass derive 
from the expectation value of a scalar field 0; the potential for depends on the number 
density of ip particles, and therefore increases naturally on cosmological timescales as the 
universe expands. As a result, the particle energy density p^ = m^n^ decreases more 
slowly than a -3 , resulting in a larger age for the universe. (There is also a contribution 
from the potential energy of 0, which redshifts at the same rate.) We discuss some of the 
cosmological consequences of this proposal, including potential observational tests. The 
question of structure formation in the presence of such particles, as well as the construction 
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of realistic particle physics models containing the necessary fields, is left for future work. 

After this paper was first submitted, we became aware of earlier an proposal for dark 
matter with time-dependent mass by Casas, Garcia-Bellido, and Quiros ||. These authors 
considered models of scalar-tensor gravity, in which the scalar coupled differently to different 
species of particles. 



2 Scale factor and age of the universe 

The model consists of a scalar and a particle species if), which can be either bosonic or 
fermionic for the purposes of this work. The mass of ip is imagined to come from the vacuum 
expectation value of 0, with the constant of proportionality some dimensionless parameter 
A: 

77ty = A(0) . (5) 

More elaborate dependences of on (0) are certainly conceivable, but for the purposes of 
this paper we make this simple choice. The dynamics of are determined by a conventional 
kinetic term and a potential energy U(<f>). The notable feature of the model is that we 
choose the potential U(4>) to blow up at = and roll monotonically to zero as —>■ oo. 
For simplicity we will write 

U{4>) = u <f)- p , (6) 

although more complicated forms are again possible. While such a potential seems unusual, 
this form can arise for example due to nonperturbative effects lifting flat directions in su- 
persymmetric gauge theories ||10H , as well as for moduli fields in string theory. (In fact this 



form of potential is not strictly necessary, as the phenomenon we will describe can occur 
with almost any potential; however, the effects are most dramatic with this choice.) 

This model possesses no stable vacuum state; in empty space tends to roll to infinity. 
We consider instead the behavior of in a homogeneous background of ?/>'s with number 
density n^. In that case, the dependence of the free energy on the value of comes both from 
the potential U(4>) and the rest energy of the i[) particles, which have a mass proportional to 
0. The equilibrium value of a homogeneous configuration is therefore one which minimizes 
an effective potential of the form 

V(<f>) = uo^'P + An^0 . (7) 
(See Fig. 1.) The additional contribution can be thought of as arising because increasing 
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Figure 1: Effective potential for 0. The light solid curve is the bare potential U(4>) oc -1 . 
The effective potential at finite density, given by the solid curves, is obtained by adding a 
contribution linear in and proportional to the number density n^. This is plotted for two 
different values of n^, corresponding to two different stages in the evolution of the universe. 
As the universe expands, decreases, and the equilibrium value of increases. 

increases the energy density in -0's since it increases the mass of ip. The expectation value 
of <t> is therefore 



(0) 



\ i/(i+p) 



(Such a configuration is not truly stable, as spatially inhomogeneous perturbations will 
tend to grow, but it can be stable enough for cosmological purposes.) Density-dependent 
potentials such as this have been discussed previously in other contexts; see e.g. |TT[ . 

In an expanding universe, the number density n$ will change with time; in turn, the 
mass of both and ip will change, as will the vacuum energy. After the interactions of i[) 
have frozen out, the number density can be written as = n^oa -3 , where n^o is the density 
when a — 1, which we take to be the present epoch. Then <fr evolves as 

= o a 3 /( 1+ ^ , (9) 

where O is the value of (H) at the present time. In terms of these variables the mass of the 
boson is given by 



ml 



d 2 V 



p(p + !H0o (P+2) a- :i ^l^ , (10) 



-3(2+p)/(l+p) 
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and the mass of ip is 

= A0 o a 3/(1+P) • (11) 

Both tp and are therefore variable-mass particles, or "vamps"; a cosmological model in 
which vamps are the dominant component of the energy density at late times will be referred 
to as VDM. 

There are a number of contributions to the energy density of the universe in this model. 
These include the energy in the scalar particles, in the ip particles, in the time derivative 
of the expectation value of 0, in the potential U ((/>), and in ordinary components of matter 
and radiation. For reasonable values of the parameters, the energies in 2 and in quanta 
are small; the former because is only changing on cosmological timescales, and the latter 
because the mass of is decreasing with time. (At early times, the bosons are very massive 
and rapidly decay.) The important new contribution is therefore simply V(<f>), the sum of 
the fundamental potential and the rest energy in the -0's. (We assume for now that ip is is 
nonrelativistic. As we discuss later, it is most likely that the particles were relativistic when 
they decoupled, but at late times their momenta have redshifted sufficiently that they are 
slowly moving today.) Both of these components turn out to depend on the scale factor in 
the same way; the ratio of the energy density in ip particles to that in the potential for is 
simply 



It is therefore convenient to deal with the sum of these two contributions, 



(12) 



(l+p\ 

p v = (l+p)u (p v = I — — I A0n</, , (13) 

which evolves as 

Pv = Pvoa- 3p/{1+p) . (14) 

The parameter w characterizing the effective equation of state of the coupled <j)/ip system is 
therefore w — — 1/(1 +p). 

The energy density pu in ordinary massive particles (baryons plus a possible cold dark 
matter component) redshifts as a -3 , more slowly than py, and will therefore be the dominant 
source of energy density in the universe for intermediate redshifts. The redshift at which 
Pv = Pm is given by 

f n \ (1+P)/3 

l + ZVM=(^) ■ (15) 

\Pmo ) 
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The age of the universe, meanwhile, will be larger than in conventional flat models. The age 
corresponding to a redshift z is given by 

<- a da' 



t 



r a da 

Jo a' 

-1 f (1+z> " r 
H / 1 — fi + ^mox + VL v $x 

Jo L 



"-+"- ,, -1/2 ( 16 ) 



(2-p)/(l+p) 



For the limiting case Qmo = 0, f2 = fVo = 1) we And that the age of the universe now is 
simply 

t = ^H, 1 (l + p- 1 ) . (17) 

Fig. 2 plots the age of universes with Q = 1 as a function of Qmo = 1 _ ^Vo, f° r P — 1- 

An interesting feature of this model, in comparison with alternative theories of rolling 
scalar fields and unusual equations of state, is that (because ip is massive and nonrelativistic 
at late times) it is at least conceivable that the energy density of the universe is dominated 
solely by baryons and vamps (without any significant cold dark matter component). For 
illustrative purposes, let us define the "minimal VDM model" as that with p — 1, il vo = 
1 — fibaryon = 0.96, and H = 70 km/sec/Mpc. This is a flat universe consisting solely of 
baryons and vamps, with the baryon density consistent with the prediction of Big Bang 
nucleosynthesis. In this minimal model, vamp-matter equality occurs at a redshift Zvm ~ 7. 
The age of the universe turns out to be approximately 17 x 10 9 years, in good accord with 
the (pre-Hipparcos) ages of the oldest globular clusters. 



3 Particle parameters and abundances 

The properties we have deduced to this point depend on the present energy density pvo, 
but not on any assumptions about the parameters of the particle physics model in which 
we imagine the necessary fields and interactions could arise. To understand the formation 
of large-scale structure in the model, however, it is necessary to know the mass and average 
velocity of the ip particles today, and to compute these requires some detailed knowledge of 
the interactions of our two fields. In the absence of a specific model, we will estimate these 
quantities under the minimal assumptions that ip was in thermal equilibrium at some high 
temperature and has evolved freely ever since. 

We begin by considering the general problem of the motion of an otherwise free particle 
whose mass = \(p may vary throughout spacetime. The motion of such a particle 
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Figure 2: Age of the universe in billions of years. The values in this plot are computed for 
flat universes consisting of only vamps and nonrelativistic matter, with p = 1. 



extremizes the action 

s = 



1/2 (18) 



where r is the proper time along the particle's trajectory and = m(dx fl / dr) is the particle's 
four-momentum. Variation of this action with respect to the path leads to an equation of 
motion 

^ = _ w ^ (19) 
dr dr p dr 

which can be written explicitly in terms of the path x^{r) as 

g? 2 x m _,, dx p dx a ( .... dx p dx v \ _ 

^ + r;-- = -^ + --)4(.„^). (20) 

Since we are assuming that <ft = <ft(t) is constant along spacelike hypersurfaces t = constant 
of the metric (]I|), we can solve explicitly for the motion of a particle obeying ([19]). In terms 
of the magnitude of the spacelike 3-momentum, 

W I 2 = QijPY = a 2 S ijP y , (21) 

we find 

W I oc a' 1 , (22) 
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just as for conventional (constant-mass) particles. The distinction arises for the velocity; if 
the four-velocity satisfies p^ = mu^ 1 , the magnitude of the three-velocity \u\ = (gijU 1 ^) 1 ^ 2 is 
proportional to (a0) _1 . Thus, as the particles get more massive with time, they naturally 
slow down even more rapidly than ordinary test particles. 

Although we have not specified any explicit interactions between the vamps and visible 
matter, we may imagine that such reactions exist as long as they are sufficiently weak that 
they do not lead to consequences which would have already been observed. As a result 
of such interactions, we presume that the ip's were in thermal equilibrium at some high 
temperature. Their equilibrium phase-space distribution function is either a Fermi-Dirac or 
Bose-Einstein distribution, 

M = | • (*» 

where is the number of spin degrees of freedom, hp is Planck's constant, k is Boltzmann's 
constant, is the temperature of the ^'s, and E is the energy, given by 

E = p° = K + \P\ 2 Y /2 ■ ( 24 ) 

As the temperature and density increase, \p\ goes up while goes down. At sufficiently 
early times, therefore, the ip particles were relativistic, and E ~ \p \ oc a -1 . Under these 
circumstances the ^'s behave like ordinary relativistic particles; their temperature redshifts 
as oc a -1 , and their energy density as oc a -4 . When they become nonrelativistic, 
on the other hand, their kinetic temperature will scale as ~ |p| 2 / m V> ^ a ~ 2 <l ) ~ 1 ', they 
cool off more rapidly than ordinary matter. (Strictly speaking it is incorrect to speak of 
a temperature after the particles become nonrelativistic, as the varying rates at which the 
particles slow down will distort the initially thermal distribution.) 

It is reasonable to assume that ip was relativistic when it decoupled, and we may proceed 
under this assumption to show that it leads to a consistent picture. In that case the number 



density of ip today is given by the standard formula | 12fl 

n^ = 825 r^p cm -3 , (25) 

where is the ratio of g e g, the effective number of degrees of freedom in ip, to the total 
effective number of relativistic degrees of freedom at freeze-out. (In terms of the number 
of physical degrees of freedom g, g e g = g for bosons and g c ^ = 3g/4 for fermions.) For 
simplicity let us consider the case p = 1. Then we can directly determine the mass of ip in 
terms of the current density parameter Q^q and Hubble constant Hq = 100/i km/sec/Mpc: 

= 12.7 tt^h 2 ^ 1 a 3/2 eV . (26) 
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In terms of the Yukawa coupling A, the other relevant parameters of the model are then 

u = 1.02 x 10~ 9 ^- (eV) 5 (27) 

and 

m 4 = 1.00 x 10- 6 -^f- a~ 9/4 eV . (28) 

The temperature of the if) particles (while they are still relativistic) is diluted somewhat with 
respect to that of the photons, due to entropy production subsequent to the freeze-out of if): 

?\ = (a*o\ 1/3 T a -i 
* WfJ 7 _ (29) 
= 3.55 x 10- 4 a~V/ /3 eV . 

Comparing (|2~6"D to (|2~9"D , we find that the ip's first become non-relativistic at a redshift of 

^ = 66.3^^^V /5 • (30) 

4 Further consequences 

Although the VDM model helps to alleviate the age problem, there are a number of other 
cosmo logical tests that could conceivably rule it out. For example, nucleosynthesis places 
stringent limits on the number of degrees of freedom contributing to the energy density at 



T ~ 1 keV |L3j . Particles which decouple at sufficiently high energies are not constrained by 
this test, as their number density is diluted by entropy production after decoupling. We do 
not know the temperature at which if> decouples, although there is no reason to believe that 
it isn't sufficiently high to evade the nucleosynthesis bound. Meanwhile, the energy density 
in cf) is much less than that in if) at high temperatures, and is therefore even less constrained. 

Any model which increases the age of the universe by changing the behavior of the scale 
factor with time will be subject to various cosmological tests which are sensitive to that 
relationship; these are conventionally used to place limits on the cosmological constant j5]. 
Currently the most promising such tests are direct measurements of deviation from the linear 



Hubble law using high- redshift Type la supernovae [14]], and volume/redshift tests provided 
by the frequency of gravitational lensing of distant quasars by intervening galaxies ||15| . 
These have recently been applied to a number of models with novel dependences of the scale 



factor on time, very similar to the scenario discussed in this paper. The results to date [16 
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seem to indicate that these tests do not rule out the kind of models considered here, but 
may be able to do so in the near future when more data is available. 

Our investigation has been exclusively in the context of an unperturbed Robertson- 
Walker cosmology The next step is to introduce perturbations and discuss CMB anisotropies 
and the formation of structure; work in this direction is in progress. However, it is worth 
noting some important features of the problem. There are two powerful effects which distin- 
guish the growth of perturbations in a VDM cosmology from conventional cold dark matter, 
and they tend to affect the power spectrum in opposite ways. The first effect is the effectively 
negative pressure of the coupled system. At zero temperature, perturbations in vamps grow 
more rapidly than those in CDM; indeed, perturbations tend to grow even in the absence of 
gravity. The other effect, meanwhile, is the free streaming of the ip particles. The ^'s decou- 
ple while relativistic, and in some respects act as hot dark matter. They will tend to flow 
out of overdense regions, damping the growth of perturbations until sufficiently late times. 
An accurate appraisal of the magnitude of this process requires numerical integration of the 
evolution equations, as the Boltzmann equation does not simplify as it would for massless 
or completely nonrelativistic particles [fT7|| . These two competing effects are not the entire 
story; for example, if the ?/>'s are fermions they will be prevented from clustering on very 
small scales by the exclusion principle |Tj| . The final perturbation spectrum is therefore the 
result of a number of processes, and cannot be reliably estimated analytically. In addition, 
of course, the simple model we have investigated here may be modified, either by altering 
the form of the potential (||) or by introducing other forms of energy in addition to baryons 
and vamps (e.g., ordinary hot or cold dark matter). 

Another direction currently under investigation is the construction of particle physics 
models in which vamps may arise. A possible origin for the scalar is as one of the moduli 
of string theory; our understanding of the nonperturbative effects which give potentials to 
such fields is not sufficiently developed to attempt realistic model building at this time. 
In supersymmetric gauge theories, however, there are (perturbatively) flat directions whose 
dynamics are somewhat better understood, and in that context the search for a model may 
be more hopeful. In such a scenario there are a number of potentially dangerous effects 
which must be avoided; for example, if the expectation value of <fi breaks supersymmetry, 
it may lead to gradual variations in the parameters of the standard model as the universe 
expands. Such variations are tightly constrained by a variety of data [|l9j . 
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